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Kinetics of pliase separation transition in boson- fermion cold atom mixtures is investigated. We 
identify the parameters at which the transition is governed by quantum nucleation mechanism, 
responsible for the formation of critical nuclei of a stable phase. We demonstrate that for low 
fermion-boson mass ratio the density dependence of quantum nucleation transition rate is experi- 
mentally observable. The crossover to macroscopic quantum tunneling regime is analyzed. Based 
on a microscopic description of interacting cold atom boson-fermion mixtures we derive an effective 
action for the critical droplet and obtain an asymptotic expression for the nucleation rate in the 
vicinity of the phase transition and near the spinodal instability of the mixed phase. We show that 
dissipation due to excitations in fermion subsystem play a dominant role close to the transition 
point. 
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I. INTRODUCTION 

Macroscopic metastable states of trapped cold atom 
systems have been a subject of active experimental and 
theoretical study for more than a decade [ij. Unlike 
a homogeneous system of bosons, where infinitesimally 
small attractive interaction between atoms leads to a 
collapse, trapped bosons are known to form long lived 
Bose-Einstein Condensates [1[3 (BEC) due to zero-point 
energy which, for sufficiently low densities, can compen- 
sate the negative interaction energy thus maintaining the 
system in equilibrium. Upon increasing the BEC den- 
sity, interaction energy grows, and, at some instability 
point (i.e., at a certain number of particles in the trap 
Nc, with Nc ~ 10'^ for a typical trap), zero-point energy 
can no longer sustain the negative pressure due to the 
interactions and the system collapses. It has been ar- 
gued in the literature 1] that near the instability point 
(for BEC densities slightly lower than the instability den- 
sity), the effective energy barrier that prevents BEC from 
collapsing becomes so low that the system can quantum 
mechanically tunnel into the dense (collapsed) state. 

Such phenomenon of Macroscopic Quantum Tunneling 
(MQT), however, has never been observed experimen- 
tally due to a strong dependence of the barrier height 
on the total number of particles in the trap (N). In- 
deed it has been shown [1] that the tunneling exponent 
for such a transition near the instability point scales as 
N(l—N/Nc)^^'^ and therefore very fine tuning of the total 
particle number N is required in order to keep the tun- 
neling exponent relatively small [(1 — N/Nc) <C 1]. Since 
for most BEC setups the total number of the trapped 
atoms fluctuates and typically obeys Poissonian statis- 
tics, the error in N scales as iV^/^ and therefore such a 
stringent requirement is hard to fulfill. Thus the system 
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is typically either in a sub-critical state with no barrier 
present (N > Nc) or is in the state with very high energy 
barrier and therefore very low MQT rate. 

In this paper we propose another paradigm for obser- 
vation of tunneling driven phase transition effects in cold 
atom systems based on the theory of quantum nucleation 
0, 0) Q- It has long been known that a mixture of '^He- 
''He undergoes a phase separation transition at relative 
concentration of ^He in ''He of around 6% at tempera- 
tures close to the absolute zero @. Since such a phase 
separation is a first order phase transition (it is observed 
to be accompanied by the latent heat release down to 
mK temperatures), the order parameter must have some 
finite (microscopic) correlation length and therefore the 
transition is expected to occur through the formation of 
nuclei of the new stable phase in the old metastable one. 
As usual, dynamics of the nucleation process is controlled 
by the competition of the surface and bulk energies of the 
nuclei and therefore, in order for a given nucleus to be- 
come stable (supercritical) , it must overcome a potential 
barrier formed by the two above contributions. While 
in most systems such a transition is a thermally acti- 
vated process, it has been argued that in the '^He-^He 
mixture at sufficiently low temperatures (below 100 mK) 
the transition is driven by the quantum tunneling. In 
particular, it was predicted 0,01 that near the transition 
line the tunneling exponent for such a transition rate is 
proportional to A/j,~^/^, where is the difference in 
the chemical potentials of the two phases. It has later 
been found experimentally that below 80 mK kinetics 
of such phase separation transition becomes independent 
on temperature and therefore it must be driven by the 
quantum tunneling. However, the experiments have been 
unable to verify the expected dependence of the nucle- 
ation rate on the systems's parameters (i.e., A/i, etc) - 
partly due to the poor knowledge of microscopic inter- 
actions between particles in such a strongly correlated 
system. 

We argue that contemporary cold atom systems pro- 
vide an excellent candidate for studying and observ- 
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ing the kinetics of such a phase separation transi- 
tion in boson- fermion mixtures @. Mixtures of bo- 
son and fermion atoms are typically realized in exper- 
iments studying fermionic superfluidity, where bosons 
play role of a coolant [9]. Another interesting realiza- 
tion of boson-fermion mixture has been demonstrated in 
two-component fermion system, where strongly bound 
Cooper pairs correspond to bosons interacting with un- 
paired fermion atoms [lo| . In the present paper we begin 
with detailed derivation and analysis of the results out- 
lined in Ref. 8. Significant attention is given to supercrit- 
ical dynamics, which reflects the dissipative mechanisms 
and is measurable in less interacting system. 



Starting from a microscopic description of a boson- 
fermion mixture we derive an effective action for the 
order parameter (the BEC density) taking into account 
fermion-boson interaction. We show explicitly that the 
classical potential for the order parameter due to such 
interaction has two minima corresponding to the two 
phases of the system (mixed and phase separated), see 
Sec, mil We analyze the coherence length associated with 
the system and demonstrate that it varies from finite to 
divergent and therefore allows different mechanisms for 
the phase transition, see Sec. IIVI At low fermion densi- 
ties the two minima of the potential are separated by the 
finite energy barrier resulting in finite coherence length, 
which points out that such a transition is of the first or- 
der [ll|, [iJI • We then derive an expression for the nucle- 
ation (tunneling) rate of the critical droplet of the pure 
fermion phase near the phase transition line and near the 
line of absolute (spinodal) instability of the mixed phase 
in Sees. IVII and IVIII We show that the transition rate is 
measurable (of the order 0.1 — Is^^/im^'^) for densities 
reasonably close to the phase transition line and small 
fermion-boson mass ratio. The limitation on the mass 
ratio comes as the result of dissipative dynamics due ex- 
citations in fermion system. Near the phase transition 
line it leads to significant modification (increase) of the 
transition rate for the system with high fermion-boson 
mass ratio, in which case the observable transition rates 
exist closer to the line of absolute instability (where the 
dissipation is less effective and the tunneling exponent 
becomes reasonably small). Our results for this regime 
are similar to those obtained in Ref. Qfor the MQT in the 
systems of trapped bosons with attractive interactions. 
However, in the case of boson-fermion separation tran- 
sition the height of the potential barrier and, thus, the 
tunneling exponent are controlled not by the total num- 
ber of particles in the trap, but their densities, scaling as 
(1 — rt/ric)^/^. Therefore, for sufficiently large numbers 
of particles in the trap the tunneling exponent can be 
fine-tuned with a desired accuracy which makes it pos- 
sible to observe the MQT rate in a well controlled and 
predictable regime. 



II. MODEL 

We consider a Bose-Einstein Condensate (BEC) in- 
teracting with a single species of fermions (in the 
same spin state). Interactions in such the mixture are 
characterized by two scattering lengths a^s and asp- 
Fermions and bosons interact through contact poten- 
tial XBpSir — r'), contributing term Xbf'4'b'4'f'^b'4'f, 
where Xbf = 2-K?i?aBF{^/'niB + l/mp); ijjB and ipF 
are boson and fermion fields respectively. In addi- 
tion, boson-boson particle interaction give rise to an- 
other term Abb?/'bV'b'/'-B''/'s/2 in the Hamiltonian den- 
sity, with \bb = 4:TTh^aBB /mB- The direct coupling 
between fermions is negligible: s-scattering channel is 
forbidden for the fermions in the same spin state, while 
p-wave scattering is small compared to boson-boson and 
boson-fermion contact interaction. The potential part of 
the energy density can be cast in the form 

Ep= - A^slV'Bp - MfV'I-V'f (2.1) 

+ ^IM" + >^bf\^b\^4'Uf 

For the purposes of the present calculation we can ne- 
glect the spatial dependence of the trapping potential and 
assume that the local densities of fermions and bosons 
are set by the constant chemical potentials iif and fiB- 
Indeed, since the nucleation occurs at finite coherence 
length (to be defined below) , the shape of the trapping 
potential should play little role in the dynamics of the 
phase transition as long as the effective size of the trap 
is much greater than the coherence length. 

It is convenient to describe the system in terms of the 
boson field only, tracing p-^/'^^^ {H is the overall Hamil- 
tonian) with respect to the fermion field. Such averag- 
ing can be easily carried out within mean field, i.e., the 
Thomas- Fermi approximation [l^, [ij, [l^, so that the 
calculation reduces to the evaluation of the canonical par- 
tition function (or free energy) of the free fermions with 
effective chemical potential jip — ^f ~ Xbf\4'b\'^ ■ Later, 
in Sec. IIVI we will proceed beyond this approximation 
and account for the fermion excitations interacting with 
the condensate. For the purpose of this and the following 
section such a correction is not necessary. 

Within Thomas-Fermi approximation [isj we deal 
with the energy density of free fermions in the zero- 
temperature limit. It can be easily found differentiating 
fermion free energy \/2Vm^J^ ji'J^ /^n'^hP' with respect to 
the volume (note that flp ^ V~^^^). We finally obtain 
the effective classical potential density for the bosons 

Eip) = ~ piBP+\xBBp' (2.2) 
- Ao (mf - XBFpf^'^ 0{pF - Xbfp), 

where e{x > 0) = 1, e{x < 0) 0, Aq = 
{2mF)^^'^ /15T:^hP , and we have used the density-phase 
variables, i.e. 4'b = y/pe^'^, for the condensate field. The 
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approximation discussed above ignores gradient terms 
in the diagrammatic expansion of the partition function 
[IB, [13] ■ This terms are not important for the phase 
structure of the equiUbrium system and will be discussed 
in Sec. |TVl 

In the next section we discuss the structure of the phase 
diagram describing possible configurations in the equi- 
librium boson-fermion system. Subsequent sections are 
devoted to the kinetics of the transition. The transition 
rates are calculated in Sees. IVII and IVIII Sections IVIIII 
and llXl are devoted to the discussion of thermal activation 
and supercritical expansion. The results are discussed in 
SecH 



III. SYSTEM IN EQUILIBRIUM: THE PHASE 
DIAGRAM 



Equation (|2.2p was obtained tracing over the fermionic 
part of the partition function. In this respect, it can be 
viewed as effective potential energy of boson condensate. 
The structure of the fast fermionic part of the system, in 
particular the density of fermions, is directly related to 
the state of the boson condensate. In the limit of zero 
temperatures we obtain the fermion density in the form 

Pf{p) = -j^^aiP'F ~ \BFpT^^d{P'F - Xbfp) (3.1) 

where the chemical potential is the same for all parts 
of the system in the mixture as far as the time scale of 
slow boson subsystem is of interest. 

It was shown within Thomas-Fermi approximation [T^ 
that the equilibrium boson-fermion system with interac- 
tions introduced above can exist in three configurations 
(phases): (i) mixture, (ii) coexistence of pure fermion 
fraction with (spatially separated) mixture, and (iii) co- 
existence of pure boson and pure fermion fractions, see 
Fig. [TJ In this section we show that the structure of the 
entire phase diagram can be obtained from the energy 
density (|2.2p . Indeed, we notice that E{p) has either 
one or two local minima. The first one (if present) is at 
p = 0. It clearly describes pure fermion fraction with 
the density of fermion pf{0). The position of the second 
minimum is found from dpE{po) ~ 0. It describes either 
the mixture with fermion density ppipo) > 0, or the pure 
boson fraction with pp[p{)) = due to the step function 
in Eq. (IXTl) . 

We first discuss the situation when the thickness of 
the transition region (boundary) between the fractions, 
or the coherence length, is microscopic (i.e. small with 
respect to the characteristic size of the fractions). In 
this case it is natural to introduce the overall (average) 
fermion pp = pf{Q)Vf + Pf{p)Vb and boson pB = VbP 
densities of a larger system. Here Vp and Vb = I — Vp 
are relative volumes for each fraction. 

When the densities ps and pp are small, the minimum 
at p = po is the only one present. At low temperatures 



the equilibrium system will occupy this minimum cre- 
ating uniform mixture (phase i) with pp = pp{po) and 
Pb = Po ("vp = 0). At higher densities the minimum at 
p — — the pure fermion fraction — forms. As a result, in 
equilibrium, the two minima align, E{0) — E{po), which 
defines the volume of the spatially separated fermion frac- 
tion. 

For clarity of further discussion it is convenient to in- 
troduce dimensionless densities [l3| as 

np = ppa%glgl and ub ^ P{B}aBB/9o- (3-2) 

Here and throughout the paper we adopt the convention 
of using "n" for dimensionless densities corresponding to 
original densities denoted by "p" . The conversion param- 
eters of Eq. p.2p are 
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(3.3) 



(3.4) 



where go is related to boson gas parameter (defined [l8[ 
via gB = \/ PBa%g) as g% = glriB- 

The separation between phases i and ii on the phase 
diagram {np vs ub) occurs when Vp sets to zero. The 
corresponding phase separation curve can be easily ob- 
tained in a parametric form as follows. Introducing 
A = Xbb /^XoX^gppy^ and = I — pqXbf/pf one can 
write the boson and fermion densities in the mixture frac- 
tion as 



y A 

"F = ^ and 



4^2 



(3.5) 



The equilibrium mixture coexists with unoccupied 
(Vp — !■ 0) pure fermion minima only at the i-ii phase 
transition curve. In terms of A and y, the equation 
E{0) = E{po) takes the form 



A 



2 + 4y + 6y'^ + 3y^ 

5(1 + vr 



(3.6) 



which, together with Eq. (|3.5p defines the phase tran- 
sition curve: n'p vs n^. From Eq. (|3.6p we also see 
that A varies within 2/5 < A < 3/4 (since by definition 
< ?/ < 1 within phase ii). 

The complete separation on pure fermion and pure bo- 
son fractions (phase iii) occurs for higher densities when 
Pa ^ Pf/Xbf- This situation is not considered in the 
present paper. 

In the discussion above, it was assumed that in the 
phase ii the fermion fraction is spatially separated from 
the mixture and the volume of the surface layer is neg- 
ligible. This is not always the case as will be demon- 
strated in the next section. However when / is large, one 
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FIG. 1: The phase diagram of the uniform boson-fermion mix- 
tures. Three phase are possible: (i) mixture, (ii) coexistence 
of pure fermions with mixed fraction, (iii) coexistence of pure 
fermion and pure boson fractions. The arrow indicates the 
path of the mixture fraction as the system equilibrates. The 
instability level is represented by the dashed horizontal line 
on the phase diagram. The densities are in units of go/a%g 
and gi/a%g for bosons and fermions respectively, as defined 
in Eq. (|3.2|l . The inset (b) shows the sketch of E{p) at equi- 
librium for the corresponding phases. The inset (c) gives the 
sketch of E{p) near (below) the instability line and (above) 
the phase separation curve. The arrow indicates (schemati- 
cally) modification of E{p) during the equilibration. 



can still define effective partial volumes Vp and Vb, as- 
suming that the densities discussed above are the peak 
densities. While the relation between Vp and Vb is now 
complicated, the limit Vp — > with Vb ^ 1 still exists. 
Therefore, the i-ii separation curve found above is valid. 



IV. DYNAMICS OF QUANTUM TRANSITION 

We are interested in the kinetics of the phase transition 
between phases i and ii, which is manifested by separation 
of pure fermion fraction out of the mixture. When the 
uniform mixture is prepared with the densities np and hb 
above the i-ii separation curve, e.g. point X, it occupies 
the second minima of E{p) and is metastable, see Fig.[T]:. 
As ni;- increases towards the point Y (and further) the 
barrier that separates local minima po becomes smaller 
and eventually disappears from np > rip^ with rip — 
8/27 (the absolute, spinodal, instability line). From this 
point the mixture becomes unstable. 

The decay of the metastable mixture results in spatial 
separation of pure fermion fraction. At low temperatures 
this transition is governed by quantum tunneling (in con- 
trast to high temperatures, when thermal activation be- 
comes effective and transition is due to thermal fluctu- 
ations rather then tunneling). To illustrate the equili- 
bration process, let us assume that the coherence length, 
I, is small enough so that the spatial separation of the 
fermion fraction is well defined. Towards the end of this 



section we will show where such condition is met. During 
the equilibration, the relative volume Vp of the fermion 
fraction increases from to some equilibrium value Vp. 
The densities po and pp{po) of the remaining mixture 
will change, as well as E{p). The point (po,Pf(po)) will 
drift towards i-ii separation curve where £'(0) = E{pq), 
as shown by the arrow in Figs. [1^ and [TJ:. As it can be 
easily verified, the drift line is parallel to ii-iii separation 
line and its intersection with the n p axis gives the equilib- 
rium value for the density in pure fermion fraction, pp{0). 
The partial volume of the later is Vp — (n^ — nB)/'n'^- 

In what follows we will obtain the density profile for the 
fermion droplet in various regions of the phase diagram 
(within phase ii) and set up equations to describe the 
transition rate in the system. The rates are calculated in 
the next sections. 

The dynamics of the boson condensate due to tunnel- 
ing part of the equilibration is given by the transition 
amplitude (M'|e-'-f^*|i') = /X'VB^'V'Be*-'^''*''""^^'^^''^^)^ 
where the rhs is Feynman's sum over the histories [T6l.lT9|. 
The state \i') represent the (non-equilibrium) metastable 
mixture residing in the second minimum of E{p), at 
p — > Po- The state correspond to the state after 
the tunneling with p — pi, E{pi) — E{po). Beyond 
this point the system hydrodynamically equilibrates to 
the state \ii) where the two minima get aligned, see 
Fig [1)3, curve ii. The Lagrangian density is defined as 
L{iPb,^*b) = ip^ihOtikB — H{iPb,iIjb)^ where the Hamil- 
tonian density contains the kinetic fi,^|VV'Bp/2mB and 
potential E(p) contributions, where the latter is con- 
sidered within the Thomas-Fermi approximation, see 
Eq. (1^ . 

As it has been pointed out earlier, approximation for 
E{p), i.e. Eq. (|2.2p . does not account for the gradient 
terms which arise due to inhomogeneity of the fermion 
subsystem and add to ?i^|V'0sP/2mB term of the La- 
grangian density. In other words, the approximation im- 
plies that renormalization of the boson kinetic energy 
arising due to the non-locality of the fermionic response 
function is relatively small. A straightforward perturba- 
tive estimate [lB| to the second order in Xbp yields the 
gradient term correction to the Thomas-Fermi of the or- 
der ~ {■m'^J'^X'^pf/J'^ /h^kp){\'p)'^. Comparing this term 
with the bare boson kinetic energy, we see that near 
the phase transition line it is smaller by the factor of 
^ ijnBl'nip^'^^^g'B^'ri^B^I'ri'p^ . Therefore, for sufficiently 
low gB and not very small fermion densities the renor- 
malization correction is negligible. This is clearly the 
case in the region near the spinodal instability line. In 
the vicinity of the phase separation line the renormaliza- 
tion correction can become important (e.g. closer to the 
tricritical point where np ^0). In this case the cor- 
rection is given by the factor of ^ {kpl)~'^, where the 
coherence length I is of the order aBB/gs or greater (as 
will be demonstrated shortly). For gB ^ 0.1 and not too 
small np (e.g. for tib ~ 0.4), (kpl)"^ ^ 50, and thus the 
Thomas-Fermi approximation is well justified. 

Considering the above arguments, the transition can 
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be described by the effective Lagrangian density of the 
form 



L — hp 



dt 



2m_B 



2mB 4/3 



E(p). (4.1) 



Here the first term can be viewed as the Berry phase; 
the second and the third terms arise due to the kinetic 
energy; the dissipative terms due to particle-whole exci- 
tations in the fermion subsystem will be included later. 
The contribution due to non-condensate component of 
bosonic system is already accounted for by the above 
treatment and therefore separate consideration is not 
necessary for the purpose of the present calculation. 

The decay rate from a metastable state can be obtained 
[TqI by calculating the classical action for the transition 
amplitude in imaginary time formalism, it — > t. Namely, 



r/F=(ro/F)exp(-5/ft) 



(4.2) 



where the action S = J dtdrL{pci, 4>ci) is evaluated over 
the classical (extremal) trajectory, Pci(0ci)- As will be 
shown below, for the parameters of interest, the rate 
is dominated by the exponent and is less sensitive to 
the changes in prefactor Tq/V. Therefore precise eval- 
uation of Tq/V is not crucial. It can be estimated as 
To /V ^ ujqI^, where / is boson coherence length as before, 
and ujQ is an "attempt" frequency. From the uncertainty 
principle loq ~ h/2mBl^ and thus Tq/V ~ h/raBl^ ■ 

The extremal action S is calculated by setting the cor- 
responding functional derivatives to zero, i.e. SL/5p = 
and 6L/S(j) = 0. The second derivative is trivial. It leads 
to the continuity equation 



dtp + V{pu)^0 



(4.3) 



with u = hS7(j)/mB and can be used to eliminate 4>. Equa- 
tion. (|4.3p can be easily solved for the velocity of the 
condensate assuming spherical symmetry. The result is 
u = (r/r^p) drr^dtp- After a straightforward algebra 
the action becomes 



S = An dtdrr"^ 



ruB f J_ 
2mB 4/9 



j drr^dtp 



(4.4) 



E{P) 



The pure fermion fraction to be formed during the 
equilibration corresponds to the bubble in the boson 
system. Its shape is defined by the last two terms of 
Eq. (14. 4p . i.e. the equation SS/5p = in the static case, 

VP + ^ a. VP - ^ d^pEip) (4.5) 

To solve it, let us first ignore the the second term of 
the left-hand side, which is a reasonable approximation 
for large enough bubbles. The solution of the remaining 
equation can be easily estimated for all relevant area of 
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FIG. 2: (a) The reciprocal coherence length as a func- 
tion of dimensionless fermion density, np, for different val- 
ues of dimensionless boson density of the mixture: ub = 
0.1, 0.2, 1.5 (from bottom to top). The dashed curve repre- 
sents reciprocal coherence length at the i-ii phase separation 
curve. 



the pase diagram, see Appendix 1X1 Near the i-ii separa- 
tion curve we obtain 



p(r) 



Po 



l + exp[-4(r-i?)/y 



where the coherence length, Z — s- is 
^0 1 



a-BB 



go^B 



(4.6) 



(4.7) 



Near the spinodal instability the first-derrivative term 
becomes important. The density p(r) varies between pi 
and pq. The coherence length I —^ Is is found, see Ap- 
pendix [XJ in the form 



Is 



V3 



a-BB 



ga^/TmB 



-1/2 



o 



The estimate of I for all the densities of phase ii be- 
low instability is given in Fig. [2^ (see Appendix |^ for 
derivations). For not too small go the characteristic 
length associated with the transition region is micro- 
scopic, i.e. I ^ CLBB- This is the case for large part 
of the interval rfip < np < rip, see Fig. [5^. There- 
fore formation of a distinct fermion fraction with thin 
boundary (nucleation) is expected for this range of pa- 
rameters. It is straightforward to obtain the character- 
istic (critical) radius, Rc of such a nuclei. Integrating 
both sides of Eq. (|4.5p with respect to y/p we obtain 
niBAE/h^ = 2 d^dr^/r. When Z < i?, we have 
Rc — a/AE, where a = {2h^ /ttib) Jq" dy/pdry/p is the 
surface tension. 



V. EFFECT OF DISSIPATION 

We now consider the effect of dissipation which was ig- 
nored so far in our description of the kinetics. Excitation 
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of the particle-hole pairs in the fermion subsystem in- 
teracting with the condensate affects the transition, as 
it happens for other systems with quantum tunneling 
El, m m HI- As will be demonstrated further, 
the resulting dissipation is important and can lead to 
significant modifications of the transition rate. 

The dissipation terms arise naturally from more accu- 
rate treatment of the boson-fermion interaction. Indeed, 
the Thomas-Fermi approximation utilized in the deriva- 
tion of the effective potential E{p), see Sees. Ulland HVl 
implies that fermions instantaneously (on a much shorter 
time scale) adjust to the local variation of boson density: 
it ignores the w-dependent terms in the diagrammatic 
expansion of the partition function. The effect can be 
analyzed by considering the second order correction pro- 
duced by the interaction \bf4'\'4'fP- 



\2 

2h? 



(27r) 



i_ -I £_!_! 

27r \47r2/^^ ^i^W-q 



(5.1) 

The first lo independent term is already present in 
Eq. (|2.2p . We are after the frequency dependent term 
rr?p\io\l ^"Kf^q responsible for Landau damping. 

In the limit I <C i?, i.e. near the phase separation 
curve and far enough from the instability line, we can use 
p(r) « po9{r — R) as suggested by Eq. (|4.6p . Substitution 
into Eq. (|5.ip gives the correction to S in the form 



AS* 



70 



dtdt' 



647r 

-R{t)R{t'f + ^[R{tf - R{t'fY In 



R{tfR{t') (5.2) 
R{t) - R{t') 



R{t) + R{t') 



Here 70 



\,p'^/Tih?. The first two terms in 



the right-hand side arise due to the restructuring of the 
fermionic density of states inside the droplet in the course 
of its expansion, while the last term can be viewed as 
coupling between droplet 's surface and particle- hole ex- 
citation in Fermi sea. 

Near spinodal instability the tunneling barrier is small 
and the correction is due to small variation of the density 
P = Po + around the metastable minima, thus 



AS- 



VI. 



■^BF 



^' ^\Sp{q,Lj)\^ 



2fi3 J (27r)3 27r Airq 



(5.3) 



NUCLEATION NEAR THE PHASE 
TRANSITION LINE 



In this section we evaluate the i-ii transition rate near 
the phase separation curve, see Fig. [T^. Evaluation of 
the extremum action with the structure of Eq. (|4.4p has 
been address in Ref. [sj. In the vicinity of the transition 
curve the coherence length, I ^ Iq, is given by Eq. (|4.7p . 
At the same time AE — > as we approach the i-ii transi- 
tion curve (the two minima of E{p) align) and, therefore, 
Rc 00. As the result Iq <C Rc and the boson density 
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FIG. 3: Dimensionless surface tension coefficient along 
the phase separation curve (numerical solution), ao = 
y/ivh^ gQ/2mBa,'BB- Approximation, Eq. (|6.3|) . is shown by 
the dashed curve. 



entering Eq. (j4.6p can be approximated by a step func- 
tion p(r, t) w pQ9[r — R{t)\ (the thin wall approximation), 
where po is bosonic density of the mixed phase as before. 
Within such an approximation the action S can be for- 
mulated [3] in terms of R{t). Evaluating the integral in 
Eq. (|4.4p over r we obtain 



S = 47r/ dt 



niBPo 



R^ {dtRf 



, (6.1) 



Here a is the surface tension defined earlier. In the 
present case (since I <^ R), the second term of Eq. (|4.5p 
is negligible and a simple energy conservation equation 
■^^(dr^/p)^ + E{p) — const corresponding to the static 
part of action S holds. Therefore we can rewrite the 
surface tension as 



2m p. 



d^p^E{p) - Eipo) (6.2) 



This integral is evaluated in Appendix |B] The numerical 
result is shown in Fig. [3) It can be approximated by 



a w 0.304 



2m f 



BO-BB 



(6.3) 



with the error of a few percents along the entire curve, 
see Appendixes E] and |B1 

It is expected that the surface tension coefficient will 
decrease as one departs from the phase separation curve. 
However near the transition line this change is not signif- 
icant compared to the change in the bulk energy. There- 
fore to estimate nucleation rates it is sufficient to use 
Eq. (|6.3p . The bulk energy is given by the integral of 
E{p) over the volume of the bubble. Therefore poAp, is 
the energy difference Appo = E{po) — E{0) or the de- 



gree of metastability, see Fig. [TJ;. When np 



A/i 



vanishes and therefore along the large part of the phase 
separation curve A/i ~ Anp. However, this in not the 



case near the end of the curve when 



0. Generally, 



Apa\ 
^BFgl 



BB 



= KAnp + K' 



3/ 2/3 01/3 



(6.4) 
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where 



K' 



2 (n 
3 



2/3 



3/2 



01/3 



K = K'^ 



2 n 



2/3 



3 ub 



-1 (6.5) 



One can see that when Anp <C nPp the second term 
in Eq. (|6.4p vanished and we can use A/j, '--^ KAnp as 
pointed out earher. In the opposite case K ^ K' ^ 
2ri^^ I'iri'p^^ and A/i ~ n^j^'^An^l''^. In this last hmit, 
however, one has to account for renormahzation of the 
boson kinetic term and, thus, the surface tension coeffi- 
cient tr, as discussed earher. 

The phase transition and formation of the fermion 
droplet is due to the interplay between the surface tension 
and bulk energy of the droplet [the last two term in 5*, see 
Eq. (|6.ip ]. which creates a potential barrier. The system 
(instanton) has to tunnel through this barrier creating 
the critical droplet of radius Rc — Sa/poAii. The rate 
of such the nucleation process is found from 63/ 6R = 0, 
with the result InT/To = -bir'^^mBPQcrRV'^ /2^/'^n'^. In 
terms of dimensionlcss densities of Fig. [TJi, for Anp <C 



Up we obtain 



Inr'/Eo = -0.0056- 



11/2 



9B 



if 7/2 \AnF 



7/2 



(6.6) 



Here is the conventional boson gas parameter as define 
earlier. 

As expected, the tunneling exponent, i.e. the rhs of 
Eq. (|6.6p . is singular in the degree of metastability Anp 

and diverges as An^^^^. Equation (|6.6p also indicates 
that the rate of nucleation is exponentially small in the 
dilute limit, i.e., for gB <C 1. Since the thin wall ap- 
proximation (nucleation) requires sufficiently high en- 
ergy barrier, e.g. Fig. [Hb,c), the rhs of Eq. (|6.6p can 
not be reduce significantly decreasing ub, see Eq. (I4.7|) 
and Fig. However, due to the smallness of the nu- 
merical coefficient one can hope that quantum nucle- 
ation is observable in sufficiently strongly coupled sys- 
tems (which are presently realizable with the use of Fes- 
hbach resonance). Indeed, for gs ^ 0.1, ub ~ 0.4, and 
Anp /rip = 0.15, the coefficient K ~ 0.27 and the tun- 
neling exponent, is ^ —27. For the same parameters and 
ubb ^ 20 a.u. the estimate of the prefactor Tq/V gives 
10^^ s~^/xm^'^. As the result the nucleation rate T/V is 
of the order 1 s~^/im~^. This is readily observable. 

Let us now include dissipation terms (|5.2p . Exact eval- 
uation of extremal action S + AS, Eqs. (|6.1|) and (|5.2p . is 
not possible and we use the variational technique. A nat- 
ural anzats is R{t) = pRce~'^^ , where coefficients a and 
p are variational parameters. The expression for AS*, 
Eq. (|5.2p . does not depend on the time scale of R[t). 
Rescaling of R{t) yields AS" ^ p'^R^. The proportion- 
ality coefficient depends only on the shape of the trial 
function, and, thus, can be estimated separately. The 
system dp{S + AS) — 0, daS = leads to the solution 



75^67^ + 42^/l07r -I- 157^1507^ + iVT^n 
3OO72 



32V157r 



(6.7) 



where for Anp ^ 



7 



0.742gl^^nB 



mp 

TUB 



4/5 



KAnp 



(6.8) 



Note that p ^^/i/^^/2 for 7 and p ^ ^/3/2 for 
7 ^ 00. The second limit gives the leading asymptotic 
in the vicinity of the phase separation curve 



InF'VFo - -0.01- 



32/5 



4/5 



mp 

TUB 



4/5 



Anp 



(6.9) 



The other limit gives expression (j6.6p with slight overesti- 
mate of the numerical prefactor due to variational nature 
of the calculation. The overall solution is 



InF/Fo _ / 8\/2 
InT'/Fo ~l TVl, 



7/2 




P+1- 



53/431/4 



-P 



(6.10) 

where we have rescaled the overall numerical coefficient 
to mach expression (j6.6p in the underdamped limit 7 — > 
0. The exponent InF/Fo as a function of Anp/np is 
shown in Fig. 2] 

The crossover between power 4 and 7/2 asymptotic 
curves take place at 7 ~ 1, or 



KAnp/ Up ^ g^ 



2/5 9/5 



mp 
ms 



8/5 



(6.11) 



Thus, we find that correction due to dissipation term, 
AS*, strongly alters the tunneling exponent in the region 

K Anp /np<^ g^^ n^J^ (mp/mB)^^^ , while in the opposite 
(non-dissipative) limit the tunneling exponent is given 
by Eq. (|6.6p . From Eqs. (|6.9p and (|6.1ip we conclude 
that the influence of dissipation is significant for large 
mp/mB mass ratio. From the crossover condition we 
also see that the dissipative regime is realized near the 
phase separation curve within a "band" of width Anp 
determined by the relation (j6.1ip . For mixtures with 
small fermion-boson mass ratio this "band" is relatively 
narrow and for large enough fermion density the system 
stays in non-dissipative regime, where the nucleation rate 
is given by Eq. (|6.6p . 

Let us consider the previous example (with gB ~ 0.1, 
nB ^ 0.4, etc.) assuming mp = mB- For this set of 
parameters the system is in the crossover region (at the 
edge of the dissipative "band" on the phase diagram). 
Using Eq. ()6.10p we obtain InF/Fo ^ —70, which means 
that the transition is not observable. Choosing a smaller 
mass ratio, e.g. mp/mB ~ 1/4 (the right-hand side of 
Eq. (j6.1ip becomes smaller by an order of magnitude), we 
shrink the dissipative region around the phase transition 
curve. The system now is in "quiet" (non-dissipative) 
regime and the previous estimates for the rate obtained 
with Eq. (|6.6p are valid. Therefore the dynamics of quan- 
tum nucleation can be systematically observable only 
for cold atom systems with sufficiently small mass ratio 
mp/mB- 
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VII. TUNNELING NEAR SPINODAL 
INSTABILITY 

In this section we evaluate the transition rate near the 
absolute instability line, i.e. rip — np^ rip. In this case 
the barrier separating the metastable mixture becomes 
small, see Fig. \Vp, curve Y. To compute the tunneling 
part of the transition based on action (|4.4p it is sufficient 
to retain only a few terms in the expansion of E{p) at 
p — > po- The terms of the second and the third order in 
5p — p ~ po will form the barrier (we set E(pq) — for 
convenience). It is clear that in the kinetic part of the ac- 
tion only the second order terms should be kept (the third 
order terms containing gradients are small compared to 
the third order potential terms). We obtain 



S = Att dtdrr 



+ 



TUB ( 1 



^(;.//.-V-av| (7.1) 



^2 i^^pf 



SrriB Po 



+ a5p^ + b6p^ 



where 



and 



a — 1 



Sgbbiti^ 



1 



mp 



The dependence of action (|7.ip on np and tib can be ob- 
tained rescaling r, i, and 5p, similarly to Ref. 3. By in- 
troducing dimensionless variables [1^ x = ry/ErriBPoa/h, 
T = Apoat/h, and p — Spb/a, the action S can be 
rewritten as const x s, where s[p{x,t)] is a parameter- 
independent functional, the extremum of which is a c- 
number. Its value can be estimated by variational anzats 
p = — Po exp (— aa;^ — /Jr^), where a, (3 and po s.re vari- 
ational parameters. We find a — /9 = 4/3, po=9, and 
s ~ 8.95. Upon a straightforward calculation one obtains 



InTy/ro 



0.324 

9Bn% 



np 



1/2 



(7.2) 



Again we see that tunneling exponent is controlled by 
the inverse boson gas parameter gB- 

The exponent vanishes when fermion density np 
reaches the instability line np, where effective energy 
barrier disappears, see Fig.[T}:. In this case the transition 
rate is defined by the prefactor Tq/V rather then expo- 
nent, which is especially troublesome for a dense system. 
This limit, however, is not of interest here. Moreover, for 
not very small nB the observable transition occurs much 
closer to the phase separation curve (far from the insta- 
bility point) and is described by Eq. (|6.10p of Sec. IVII 

When gB S 0-01 and ns ~ 0.1 — 1, the phase transition 
is not observable near the transition curve, see Eqs. (j6.6p 
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FIG. 4: Transition rate exponent, InF/Fo, as a function of 
fermion density, np = n% + AnF, for tib = 0.4, qb = 0.1 (left) 
and gB ~ 0.01 (riglit). Asymptotics near the i-ii phase sep- 
aration curve are plotted for mp/mB ~ 0.01, 0.1, 1, 10 (solid 
curves, bottom to top). The dashed line represents the rate 
without dissipation, see Eq. (|6.6|l . The asymptotic near the 
instability region, Eq. (|7.2|l . is given in the right plot (for the 
same variation of mp/mB)- 



and (|6.9p . — one has to chose the values of n^^- closer to 
the spinodal instability. In this case the exponent (|7.2p 
is significant. It rises rapidly as the function of np ^ np 
defining the transition rate near the instability region. In 
this respect Eq. (|7.2p is similar to the results on Macro- 
scopic Quantum Tunneling (MQT) in systems of trapped 
bosons with attractive interactions ,lj. When the coher- 
ence length Is, see Eq. (|4.8p is still much smaller then the 
size of the condensate, the height of the potential barrier 
and thus the tunneling exponent for the critical droplet 
are controlled not by the total number of particles in the 
trap, but the local densities. Therefore, for sufficiently 
large numbers of particles in the trap the tunneling ex- 
ponent can be fine-tuned with a desired accuracy, which 
can make observation of the MQT rate possible in a well 
controlled and predictable regime. 

The contribution of dissipation, e.g., action AS* in 
Eq. (|5.3p , can be estimated by using the same variational 
anzats as above. Rewriting the integral in Eq. (j5.3p in 
terms of x, t and p one finds that AS* does not depend on 
a and thus np — np. Therefore for large enough np — rip 
and small gB one can treat the dissipation perturbatively 
substituting p{x, t) from the above calculations. As the 
result the tunneling exponent acquires an additional term 

"(S-B't^b't^f/'tib)^^^- This term is again controlled by 
the mass ratio mp/mB. For not too high fermion/boson 
mass ratio, it is of the order gB and thus dissipation does 
not significantly alter the dynamics of the phase transi- 
tion (MQT) in this region, see Fig. S) 



VIII. THERMAL NUCLEATION 

It is instructive to analyze the influence of the classi- 
cal (thermal) activation mechanism. In order leave the 
metastable state finite temperature system can (in addi- 
tion to direct tunneling leakage from its false "ground" 



9 



state) use its entire spectrum tunneling from high en- 
ergy states or going over the barrier. The rate of 
such transition is proportional to the sum [3, [2^ of 
r(£'„) exp[—{En — Eo)/kT], where the prefactor accounts 
for possible tunneling from n-th exited state of the 
metastable mixture. This regime becomes effective when 
the largest energy gap — between the first exited and the 
ground state energies of the mixture — is of the order 
ksT. This energy difference enters the expression for 
quantum nucleation rate and can be easily estimated in 
a similar manner 



El — Eq 



P 



1 



2mB 2mB \ I 



Here the characteristic length scale is the coherence 
length of boson-fermion mixture. 

Near the phase transition curve I is given by Eq. (|4.7p . 
i.e. I ^ a-BB / 9on-B- This gives the relation 

glriB - I2 (8-2) 

Critical temperature of BEC transition can be estimate 
by that of the uniform non-interacting three-dimensional 
boson gas fesTc ^ h^n'g^ /ruB, and we finally obtain con- 
dition for thermal assisted transition in the form 



2/3 J 

9b iT'B'^ y 



(8.3) 



Our boson system is in deep BEC state, i.e. T ^ Tc, 
therefore thermal fluctuations are not significant for ob- 
servable quantum nucleation that occurs in sufficiently 
dense systems {ub ~ 1). 

In dilute systems, where direct tunneling from low lay- 
ing states is extremely slow, one has two compute the 
entire sum 



(8.4) 



up to the continuum states where 5„ — > 0. The sum- 
mation can be carried out semiclassically minimizing the 
new energy functional, i.e. exponent of Eq. (|8.4[) . see e.g. 
Ref. i 

Closer to the spinodal instability the barrier becomes 
smaller and the tunneling is effective even for relatively 
dilute systems. For the crossover temperature we obtain 



2/3 



(1 - np/n%) ~ T/T„ 



(8.5) 



In this expected, the zero-temperature MQT 

mechanism can be observed in certain distance from the 
instability line, i.e. when the left-hand side of Eq. 
is greater. 



IX. EXPANSION OF SUPERCRITICAL 
DROPLET 

As soon as a critical droplet is created as the result 
of tunneling or temperature fluctuations, it should grow 



equilibrating the system as pointed out in Sec. lIVI Below 
we discuss only the case of finite coherent length — the sit- 
uation when a distinct stable droplet is created (thin-wall 
approximation). We will estimate the expansion rate for 
I ^ Rc R in the limits of large and small dissipation. 

After a given nuclei passes (tunnel trough) the energy 
barrier, its growth can be described by classical equa- 
tion of motion for radius R, 6S/6R + fd = 0- Here fd 
is the dissipative force (external time-dependent force is 
negligible at low temperatures), and the action S should 
be formulated in real (original) time, i.e. in action (j6.ip 
one replaces it —^ t. In the underdamped limit, fd^Q 
(energy conserved), the expansion is described by 



ruBPoR^ 



{dtRY 



-(jR^ + A/i 



PoR^ 



(9.1) 



For large droplets {R ^ Rc) the surface tension term is 
negligible and the expansion is linear in time R ~ agt, 
with 



y^2Ap/3mB- 



(9.2) 



Fermion particle-hole excitations affect the the expan- 
sion dynamics, renormalizing the expansion rate ag- To 
estimate this effect we need to account for the dissipa- 
tive part of the action. This can be done by analytical 
continuation of Matsubara action, Eq. (|5.ip . to Keldysh 
contour. Following Ref. [13 we obtain 



m2 \2 



'^i^^^^Pc/('7, t)pq{q, t) + 0(p^(g, t)). 



(9.3) 

where Pd/q = {p+ ± P-)l "/^ and p± resides on forward 
(backward) branch. In the nucleation limit (thin wall 
approximation) the time dependance of p(r, t) enters only 
via radius of the droplet, R{t), therefore p± = p{R±\. In 
the classical limit we keep only terms linear in Rq. The 
dissipative force becomes 



/ = 



dt 



dq 

(2^3 q 



1 • / .JPq{<ht) 

-pcM,t)-^ 



5R 



(9.4) 



Substitution of p{r,t) — po6{r — R{t)) will lead to the 
logarithmic divergence in the integration over the mo- 
mentum. The divergence is due to the high q. The cut- 
off comes from the thickness of the surface area of the 
droplet. Using a more accurate expression for p{r,t), i.e. 
Eq. (j4.6p . we obtain 



/ = -jeloR{t)Rity In 



2, 4i?(t) 



ttI 



(9.5) 



where 7e — 0.577... is the Eulcr constant, and I is the 
coherence length. As the result, expansion of a large 
supercritical droplet is governed by 



niBPoRdfR 



SruBPo 



[dtR? 



(9.6) 



+ 7e 70 (9ti?) In 



4i? 
Id 



poAp. = 0. 
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This equation can be analyzed numerically for dimen- 
sionless variables. Define t — tQT, R{t) — Rcx{t), then 
Eq. (|9.6p becomes 

- Xdh - lidrxf - UdrX) ln(x/0 + 1 = 0. (9.7) 

Here /o = jeloto/mBPoRc, £. = ttI/4Rc, and = 
R^tub /A/1. Numerical integration shows that for /o ^ 1 
and T > lOO/o the first and the second terms remain 
smaller by several orders of magnitude. Therefore we 
can ignore this terms in such the limit. Moreover, we 
notice that logarithmic pre-factor of the third term does 
not change significantly for large supercritical droplets. 
Therefore with sufficient accuracy the growth is still lin- 
ear in time R{t) ~ athi~^ {Aat/irl) with 



7e70 

One can notice that fo ^ l/y/Ajj, and therefore di- 
verges at the phase separation curve. As the result, the 
above solution is valid near the phase separation line. 
When the system is placed farther away from the phase 
separation curve, A/i increases and /o becomes smaller. 
In this case we have to retain the first two derivative 
terms in Eq. (|9.7p . It changes the asymptotic behavior. 
In the limit /g ^ we arrive to Eq. (|9.ip with the solu- 
tion R{t) aot, where ao is given by Eq. (|9.2p . 

Similar to nucleation, the expansion of a supercriti- 
cal droplet is controlled by particle-whole excitations if 
the system is prepared within a certain "dissipative" re- 
gion around the phase separation line, e.g., Sec. IVII For 
larger degree of supersaturation, dissipation is less ef- 
fective. The "dissipative" region can be approximately 
estimated from condition /o(A/i) > 1. We should note 
that measurement of the growth rate does not necessar- 
ily require extremely dense systems in which nucleation 
is observable. Supercritical droplet can be created ex- 
ternally, e.g. by some laser fields or trapped atoms of 
different sort. For such initiated supercritical droplet the 
growth rate could be readily observable. 

X. DISCUSSION 

In the previous sections we have considered dynamics 
of the nuclei of a new phase (i.e., pure fermion phase) in 
the fermion-boson mixture in the metastable state. We 
have derived asymptotic expressions for the nucleation 
rates for the system in two regimes: in the regime of 
weak metastability, e.g., near the phase transition line 
(Sec. IVI|) and near the spinodal (absolute) instability. 

Analysis of Sec. IVII leads to the conclusion that in the 
vicinity of the phase transition line the nucleation dy- 
namics is practically not observable for sufficiently di- 
lute systems, e. g., with gB <C 1. It may appear 
that the dependence of Eqs. (|6.6p . (|6.9p . and (|6.10p on 
Ub contradicts to this statement. Indeed the transition 



rate decreases with ns if Anp is kept constant. Such 
limit, however, is incorrect, since the nucleation condi- 
tion ^0 Rc in not satisfied, Iq ^ I/ub (or \/ ^ub if 
qb is constant). Evaluating Rc in terms of Anp one 
obtains the condition on Anp for which the asymptotes 
(HH), dnH), and (PrU)) are valid. We have Rc/ubb = 

3 /2 

0. 219{ng /gBK){np/ Anp)- In the limit ^ it gives 
Anp/rip ^ nB (where the densities are dimensionless, 

1. e. in units of Fig. [T^, as before). Therefore in ub — > 
limit the power 7/2 and 4 asymptotes are not measurable. 

Measurable transition rates for dilute systems can be 
found closer to the instability region. In this case one 
enters a crossover between nucleation and MQT regimes. 
As an example, let us consider the mixture with the bo- 
son gas parameter qb ~ 0.01. For ubb ~ lOnm and 
about 10^ confined particles. This corresponds to the 
densities of the order lO^^'cm"'^ with size of the trapped 
condensate ^ lOO/im. For this parameters the posi- 
tion of the metastable mixture on the phase diagram. 
Fig. [T^, is primarily defined by the ratio aBs/ciBF, see 
Eqs. p.3p and (|3.4p . and can be controlled by chang- 
ing the scattering lengths. The transition rate is T/V ^ 

exp(— x/gsn^) where x"^ = l — np/up, see Eqs. (|7.2p . 
(|4.8p . The change of the prefactor is not significant due 
to fast decay of the exponent form x > O.ln^ (for smaller 
X accurate evaluation of Tq/V is necessary). Therefore 
we conclude that for {np — np)/np > O.Oln^ and ub ^ 1 
the dependence of the exponent (|7.2p on the densities is 
measurable. In this case the coherence length Is ^ l/im. 
It rises for small rip — np, however at Is <SC 100/im (or 
1 — np/up ^ 10^"*) the transition is still controlled by 
local densities. 

We have also found, e.g. Sections IVl and IVll that dis- 
sipation can rather strongly alter the dependence of the 
nucleation rate in the vicinity of the phase transition line. 
While the results obtained in these sections are likely 
to be correct only quantitatively, they show that in the 
vicinity of the transition line excitation of particle-hole 
pairs in the Fermi sea plays crucial role. As a result 
the Thomas-Fermi approximation is inapplicable within 
the "band" controlled by the fermion-boson interaction 
strength as well as fermion-boson mass ratio; see Sec. IVII 

Finally we considered the dynamics of supercritical 
droplets, i.e, the ones which have tunneled to the super- 
critical size. We have found that in the non-dissipative 
regime the radius of the droplet grows linearly with 
time with the expansion rate coefficient proportional to 
y/Anp. In dissipative region, which occurs near the 
phase transition line, the situation is quite similar: up 
to logarithmic corrections the radii of the droplets still 
grow linearly with time, with rate now being dependent 
on friction coefficient, see Sec. IIXI 
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APPENDIX A: DENSITY PROFILE 

To analyze the shape of the condensate density func- 
tion p(r) of the bubble, i.e. the solution of Eq. (|4.5p . 
it is convenient to introduce dimensionless quantities as 
X = r/r, ^ = t/t, £ = E/<B, s = S/&, u = pXbf/^J-f, and 
write the action (14.411 in the form 



/dS^x^dx — ( — 77 / x^dfudx 
u \ x^ J 



(dxu) 



elu 



with 



e{u) 



5Auo + -il 



3/2 



5 . 2 

2 



(1 



(A.l) 



N5/2 



(A-2) 

where uq = XbfPo/ij-f = 1 - J/^, t = H/AXqXbfIJ'f^, 
= ht/2'mB, and € = Xq^'^J^ . Note that the parameters 
A and y (or uo), introduced in Sec. IIIH are generally 
independent and define the position on the phase diagram 
via the same relation as in Eq. p.5[) . The relation (|3.6[) 
between them defines the special case — the equilibrium 
state (residing on the phase separation curve). 

Here we are after the time-independent part of the ac- 
tion s. In the current notation equation ()4.5p reads 



1 de{u) 



4 dy^ 



(A.3) 



Let us first discuss the case when the second term of the 
left-hand side is negligible (which is clearly the case near 
the phase transition curve). Equation (jA.3|) can, then, 
be simplified to 



(A.4) 



noticing that {dx\/u)d^ = d^- This is also the conse- 
quence of energy conservation at large x corresponding to 
the static part of action s. The right-hand side of (|A.4p is 
too complicated for the equation to be solved analytically. 
However, for the purpose of integration, e{u) — e(uo) 
can be approximated in the range ui < u < uo, where 
e{ui) — e{uo) — 0, as 



^Je{u) - e{uo) « m{u - ui)^/^(uo - u) (A.5) 



with 



maximum of the right-hand side of Eq. (jA.4[) and there- 
fore should lead to correct estimate of the characteristic 
length scale of u{x). The maximal absolute deviation 
I7(m) = l — 4^$H^(u — ui)(uo — u)7 [£(«) — e(Mo)] occurs at 
M ^ Mo and u — > Ml at the end of the i-ii separation curve 
(mq — > 1), as can be easily verified numerically. At that 
point we obtain t]{uq) « 0.421. The error vanishes at the 
limit of small barrier, i.e. near the instability line on the 
phase diagram. With this approximation Eq. (IA.4|) can 
be easily solved with the result 

exp \AW^Juo\/uQ — ui{x — const)\ = (A. 7) 

_ {^/u + y/u^) {y/UUQ + yju - Uly^UQ - Ul - Ul) 



UlV^O - Ml - Ml) 



Here the characteristic size of the bubble, R, enters into 
the constant of integration. The characteristic length- 
scale on the surface of the bubble (the coherence length) 



(A., 



2Umoa/1 - mi/mo 



where v/ubb = ^JhflAc'l'^ j \/2jig^. 

In the vicinity of the i-ii separation curve Ui — > 0, 
A is a function of mo according to Ea. (|3.6p . and we use 
Urn — uq/2. Substitution to Eq. (jA.6p yields 



5/2 



2m3 



15 A 

16 Mo 



5(l-Mo)3/2 



(A.9) 



As the function of mq, WVA varies from 0.1195 to 0.1212 
with the mean value of 0.1203 (the root-mean-square de- 
viation is 0.0004). Therefore near the i-ii separation curve 
with the same precision the coherence length, Iq, becomes 



^0 



1 



go^B 



(A.IO) 



In the vicinity of spinodal instability uq—ui ~ 2^ A^l\— 
riF/np)/^ verified by expanding e in powers of 

Mq — Ml, and we use Um — (mo + 2mi)/3. Substitution into 
Eq. jO)) yields 



2U^ 



5mo 



256V1 - Mo 



+ C'(mo - Ml) 



(A.ll) 



Therefore the coherence length, I 



obtain 



clearly diverges as 
riF/np- Collecting all dimensional factors we 



x/3 



O-BB 



gOy/TmB 



^ riF 



-1/2 



-O 1 



2^2 ^ £(""0 - ^("o) 



42 (i 



Mi)(mo - M„i)2 



Here mi < m,„ < mq and m,„ ~ (mi 4 
approximation is exact at m = Um, i.e. 



(A.6) 

Mo)/2. This 
close to the 



In the above analysis we have ignored the first deriva- 
tive term in Eq. (|A.3p . It is strictly speaking negligi- 
ble only near the phase transition curve for large thin- 
boundary bubbles. In the instability region solution 
(lA.yp will change. Nevertheless, the analysis of Eq. (|A.3p 
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2 0.6 1 

z 



where we used the notations of Appendix [X] Utihzing 
the approximation (jA.Sp at the i-ii transiiton curve, i.e. 
setting u\ = 0, we obtain 



d^^e(uQz) + 1 = 22Jwo/^ [1 + ©(ns/lOO)] (B.2) 



FIG. 5: Normalized energy deviation, see Eq. (IB.2[) . for uo ~ 
0.1, 0.2, 0.9. Inset (a) gives absolute error of this expression 
with respect to a straight line approximation (for the same set 
of uo). Inset (b) shows the relative error between the result 
of numerical integration in the Ihs of Eq. HB.2P and SZHiig''^ as 
a function of ns (solid line). The dashed line is a linear fit of 
the form O.OSSns. 



The left-hand side can be integrated numerically for var- 
ious uo to verify the result, see Fig. [51 Significant relative 
error of Eq. (jA.Sp at u ^ uq is suppressed in the above 
integral and the approximation (|B.2p is accurate along 
the entire phase separation curve. Using the value of 371 
found earlier (see Appendix |X| we finally obtain 



shows taht the characteristic length scale still has to be 
inverse-propo rtional to Mmq and incorporate the diver- 
gence of l/\/l — ui/uq. Therefore relation (jA.8|) still r^^l^ 

holds. This can be demonstrated by expanding the right- cr = 0.304— — [1 -t- 0(1/100)] (B.3) 

hand side of Eq. (|X3|) in terms of 1 - mi/uo. In the 2mBagg 
intermediate region, i.e. < uq ~ ^^i < ^^Oj the order-of- 
magnitude agrement is expected since no other features 
are present, as can also be verified numerically. 



APPENDIX B: SURFACE TENSION NEAR THE 
TRANSITION CURVE 

The surface tension has been defined as 
* 



For large ub, closer to "tricritical" point, the effective 
surface tension will deviate from the one obtained above 
due to additional gradient terms in Eq. (|4.4p . This is 
due to renormalization of the boson kinetic energy which 
is ignored by the Thomas-Fermi approximation. This 
correction is negligible for the parameters of interest, as 
explained in Sec. II VI 
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